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^INFLUENCE  DE  LA  CINETIQUE  CHIMIQUE  SUR  LA  COMBUSTION 

DES  GOUTTES 

fttsumt 

L'influence  de  la  cinltique  chimique  sur  la  combustion  des  gouttes  est  r&lisle  avec 
un  module  qui  considire  sym6trie  sph6rique  et  les  conditions  quasi-stationnaires.  La 
cin&ique  chimique  du  prods  est  abordd  avec  une  reaction  totale  d’ordre  n. 

Les  resultats  son  exprim&  en  fonction  du  produit  pn/2  r„,  p  6tant  la  pression, 
r„  la  rayon  de  la  goutte  et  n  1’ordre  de  la  reaction  totale.  On  dlduit  qu’il  existe  une 
valeur  minimum  de  ce  produit  en  dessous  duquel  une  goutte  isolfc  ne  peut  pas 
maintenir  une  flamme. 

On  a  effectue  aussi  une  recherche  explrimentale.  On  a  observe  la  combustion  dans 
l'air  i  pression  variable  de  gouttes  de  combustibles  suspenducs  par  des  fils  de  quartz 
tres  minces.  La  variation  avec  le  temps  du  diam&tre  de  la  flamme  et  de  la  goutte 
furent  obtenus  ainsi  que  les  valeurs  minimums  des  diamitre  des  gouttes  pour  bruler. 

On  montre  que  les  r&ultats  explrimentaux  pr6sentent  un  accord  qualitatif  excellent 
avec  ceux  de  la  thlorie. 

Abstract 

The  influence  of  chemical  kinetics  on  the  combustion  of  fuel  droplets  is  studied 
by  means  of  a  theoretical  model  of  the  process  which  considers  spherical  symmetry 
and  quasi-stationary  conditions.  Chemical  kinetics  is  approximated  by  means  of  an 
overall  reaction  rate  of  n01  order.  r_ 

Results  are  expressed  as  a  function_of  the, prod uctfp^^  r,,  where  p  is  the  pressure,  t  r  , 
n  the  order  of  the  overall  reaction  and  r^the  droplet  radius.  It  is  shown  that  there 
exists  a  minimum  value  of  such  product  under  which  an  individual  droplet  cannot 
sustain  a  flame. 

An  experimental  investigation  was  also  carried  out,  by  burning  in  air  at  variable 
pressure  fuel  droplets  suspended  of  thin  quarts  fibres.  The  laws  of  variation  of  both 
flame  and  droplet  diameter  as  functions  of  time  were  obtained  as  well  as  the 
minimum  values  of  the  droplet  diameters  for  combustion  as  a  function  of  pressure. 


/ 
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It  is  shown  that  all  experimental  results  are  in  excellent  qualitative  agreement  with 
those  predicted  by  theory.  (.  . _ 


1.  INTRODUCTION 

The  complex  combustion  phenomena  which  take  place  within  the  combustion 
chambers  of  gas  turbines  or  compression-ignition  engines  are  essentially  based  on 
the  process  of  the  combustion  of  a  spray  of  different  sized  droplets,  which  are  in 
motion  within  an  oxidizing  atmosphere. 

Combustion  of  sprays  is  studied  by  means  of  statistical  methods,  by  applying  the 
laws  of  the  combustion  process  of  individual  droplets  in  a  still  and  oxidizing  atmos¬ 
phere.  This  process  is  modified,  but  not  essentially,  i>?  by  the  proximity  of  other 
burning  droplets  (interaction  effects)  and  by  the  motion  of  the  droplets  within  the 
oxidizing  atmosphere  (forced  convection  effects). 

Combustion  of  individual  droplets  has  been  generally  studied  by  assuming  that 
the  chemical  reaction  rate  is  infinitely  fast.  This  is  the  Burke  and  Schumann  hy¬ 
pothesis  for  studying  diffusion  flames,  which  implies  that  chemical  reaction  occurs 
at  a  zero-thickness  layer,  where  the  mass  fractions  of  both  fuel  and  oxidizer  become 
equal  to  zero.  Chemical  kinetics  is  disregarded  and  results  are,  then,  easily  obtained, 
but  it  will  be  shown  that  this  assumption  may  introduce  important  errors  in  the 
analysis  of  the  process. 

On  the  basis  of  the  aforementioned  assumption,  it  is  obtained  that  the  burning 
rate  of  the  droplets  is  proportional  to  the  droplet  radius  and,  therefore,  that  the 
square  of  the  droplet  radius  is  a  linear  function  of  time.  It  is  also  obtained  that  the 
flame  temperature  and  the  flame/droplet  radius  ratio  do  not  depend  on  the  droplet 
radius.  Another  important  conclusion  is  that  pressure  does  not  influence  the  process. 

There  is  available  a  large  experimental  evidence  on  the  combustion  of  droplets. 
Experimental  techniques  consisted  in  observing  the  combustion  of  droplets  sus¬ 
pended  of  thin  quartz  fibres  *>2,  or  simulating  the  combustion  of  liquid  droplets  by 
means  of  porous  spheres  3.4  to  which  a  controllable  fuel  flow  is  supplied.  There 
also  exists  a  few  experimental  data  obtained  with  ejected  and  free  falling  drop¬ 
lets  5,0. 

Droplets  of  practical  interest  are  smaller  than  200  n  in  diameter,  as  shown  in 
Fig.  1.  It  is  very  difficult  to  observe  the  individual  combustion  of  such  small 
droplets.  Therefore,  the  largest  part  of  the  available  experimental  evidence  has 
been  obtained  with  droplets  larger  than  0.;  mm  in  diameter.  Experimental  results 
obtained  with  large  droplets  have  been  extrapolated  for  small  ones,  but  it  will  be 
shown  in  this  work  that  such  extrapolation  is  not  valid. 

The  experimental  values  of  the  burning  rates  of  large  droplets,  measured  by  the 
2 

slope  of  curves  r  =  f  (t),  agreed  fairly  well  with  their  theoretical  values,  although 

Hall  and  Diederichsen  and  Bolt  5.6  found  experimental  values  for  small  droplets 
smaller  than  those  obtained  for  large  droplets. 

The  theoretical  and  experimental  values  of  the  flame/droplets  radius  ratio  did 
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not  coincide,  and  the  experimental  values  were  not  constant,  but  increasing  with  the 
droplet  radius.  Such  disagreement  was  attributed  to  the  influence  of  free  convection, 
but  Kumagai  and  Isoda  found  similar  laws  of  variation  in  a  zero-gravity  field  7>8. 
They  explained  the  difference  between  theoretical  and  experimental  results  by 
means  of  an  approximate  theory  on  transient  combustion  of  droplets.  However, 
transient  phenomena  alone  cannot  produce  such  effects,  because  Wise,  Lorell  and 
Wood  3  found  similar  laws  of  variation  of  the  flame/droplet  radius  ratio  for 
steady-state  combustion  by  using  porous  spheres  of  different  sizes. 

It  will  be  shown  that  such  disagreements  between  theoretical  and  experimental 
results  are  due  to  the  theoretical  model  of  the  process,  because  the  assumption  of 
considering  an  infinitely  fast  reaction  rate  is  only  strictly  applicable  when  either  the 
droplet  radius  or  die  pressure  tend  towards  an  infinite  value. 

There  are  very  few  theoretical  studies  on  the  combustion  of  droplets  considering 
finite  chemical  kinetics.  Lorell,  Wise  and  Carr  ®  studied  the  problem  by  taking  an 
overall  reaction  rate  and  solving  numerically  the  differential  equations  of  the 
process.  However,  they  took  a  constant  and  rather  large  value  of  the  droplet  radius, 
they  did  not  consider  the  influence  of  pressure  and  they  studied  the  problem  for 
the  case  of  equal  molecular  weights  of  all  chemical  species.  Under  such  conditions 
results  obtained  by  assuming  an  infinite  reaction  rate  are  not  very  different  from 
those  derived  considering  finite  chemical  kinetics. 

At  the  Institute  Nacional  de  T&nica  Aeroniutica,  Madrid,  Spain,  a  research 
program  on  combustion  has  been  conducted,  part  of  which  comprised  the  study  of 
the  influence  of  chemical  kinetics  on  the  combustion  of  liquids.  In  the  first  place, 
a  general  theoretical  study  was  carried  out  on  the  combustion  of  droplets  within  a 
reactant  atmosphere  considering  second  order  chemical  kinetics  10.  Afterwards,  the 
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theory  was  applied  to  the  case  of  the  combustion  of  bromine  droplets  in  a  hydrogen 
atmosphere  u,  considering  the  actual  chemical  kinetics  of  the  process. 

At  the  same  time  an  experimental  investigation  was  performed.  The  investigation 
consisted  in  observing  the  combustion  of  bromine  and  nitric  acid  droplets  in 
hydrogen  u  at  variable  pressure  and  in  observing  the  combustion  of  several  types 
of  fuel  droplets  in  air  at  low  pressure. 

This  investigation  was  performed  in  order  to  verify  the  theoretical  conclusions 
and  in  order  to  obtain  some  practical  data  which  might  be  applied  to  the  study  of 
the  combustion  of  fuel  sprays. 

In  the  present  work  a  brief  r£sum£  of  the  theoretical  study  of  the  combustion  of 
droplets  considering  finite  chemical  kinetics  will  be  included,  as  well  as  the  prin¬ 
cipal  results  which  have  been  obtained.  Experimental  results  obtained  by  burning 
fuel  droplets  in  air  at  variable  pressure  will  also  be  included,  as  well  as  a  qualitative 
comparison  between  theoretical  and  experimental  results. 

It  will  be  shown  that  chemical  kinetics  influences  decisively  the  process  at  low 
pressure  or  for  small  droplets  and  that  the  assumption  of  considering  an  infinitely 
fast  reaction  rate  is  an  asymptotic  solution  of  the  process,  which  is  the  limiting  case 
of  the  actual  solution  of  the  problem  when  either  the  pressure  or  the  droplet  radius 
tend  towards  infinity. 

An  important  conclusion  of  the  theoretical  study  considering  finite  chemical 
kinetics  is  that  there  exist  minimum  values  of  either  the  pressure  or  the  droplet 
radius  under  which  an  individual  droplet  cannot  sustain  a  flame.  This  interesting 
conclusion  has  been  experimentally  confirmed. 

2.  FUNDAMENTAL  ASSUMPTIONS  AND  GENERAL  EQUATIONS*) 

The  theoretical  treatment  of  the  process  is  based  on  the  following  assumptions: 

The  droplets  are  considered  to  be  isolated  and  at  rest,  and  the  process  is  assumed 
to  have  spherical  symmetry.  Therefore,  the  influence  of  natural  convection  is 
disregarded.  The  process  is  stationary.  Under  such  conditions  the  study  is  only 
strictly  applicable  to  the  combustion  of  constant  radius  droplets  fed  with  a  fuel 
flow  equal  to  the  amount  of  fuel  which  is  evaporated  and  burned. 

However,  it  has  been  shown  that  for  large  or  medium  sized  droplets  the  errors 
introduced  by  considering  steady-state  conditions  are  not  important.  This  is  due  to 
the  fact  that  radial  velocity  of  the  droplet  surface  is  small  as  compared  to  the 
diffusion  velocities  of  the  species.  For  this  reason,  results  obtained  by  using  porous 
spheres  of  constant  diameters  are  similar  to  those  obtained  by  burning  real  liquid 
droplets. 

In  order  to  obtain  general  conclusions,  only  reactant  species  and  reaction  products 
will  be  considered,  and  the  actual  chemical  kinetics  of  the  process  will  be  approxi¬ 
mated  by  means  of  an  over-all  reaction  of  second  order.  However,  the  analytical 

*)  This  r6sum£  of  the  theoretical  treatment  of  the  process  has  been  taken  from  references 

10  and  12. 
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method  used  to  solve  the  problem  is  applicable  to  more  complicated  kinetic 
schemes. 

For  simplicity,  average  values  will  be  adopted  for  the  thermal  conductivity  and 
for  the  specific  heat  of  the  mixture.  The  gas  pressure  will  be  considered  constant 
throughout  the  process  and  the  heat  transferred  to  the  droplet  surface  through 
radiation  will  be  neglected. 

Under  the  aforementioned  assumptions  the  general  equations  of  the  process  are 
as  follows : 


a.  Continuity 

If  w,  is  the  reaction  rate,  fj  the  flux  fraction  of  species  i,  m  the  droplet  burning 
rate  and  r  the  radius,  we  have  for  each  chemical  species : 


m  d  f, 
4  n  r2  dr 


0) 


If  fI{  and  »'p  are  the  stoichiometric  coefficients  of  the  reactants  and  products  of  the 
over-all  reaction,  and  MR  and  MP  the  molecular  masses,  we  have : 


or  else: 


WR  __ 

wP 

vr  Mr 

yp  Mp 

it  is  obtained: 

fR  “  CR» _ 

*p  —  £p* 

’’r  Mr 

Vp  Mp 

f Roo  CR»  __ 

_fPoo—  *P* 

VR  Mr 

pp  Mp 

(2) 

(3) 

(4) 


b.  Energy 

If  hj  is  the  specific  enthalpy  of  species  i,  T  the  temperature  and  A  the  average  value 
for  the  mixture  of  the  thermal  conductivity,  the  equation  of  energy  is  given  by : 


m  2  h,  c,  —  4  a  r2  X  -j-  =  cte  (5) 

which  expresses  that  the  sum  of  the  fluxes  of  enthalpy  and  heat  is  constant.  The 
value  of  this  constant  is  obtained  by  expressing  that  at  the  droplet  surface  the  heat 
transmitted  is  used  in  vaporizing  the  liquid  fuel,  that  is : 

*  dT 

4  *  f2.  ("  jjr).  =  (6) 

where  q,  is  the  latent  heat  of  evaporation  of  the  fuel. 

It  results : 
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m 


^hi.*i.  +  qi)-4«r«T^E  =  o 

i  i  dr 

Referring  the  enthalpies  h,  to  the  standard  state  hj°  in  the  form : 

h,  =  V  +  <„«  (T  —  T0),  (8) 

taking  an  average  value  cp  of  the  specific  heats  for  the  mixture  and  introducing  the 
heat  of  reaction  qr  by  means  of  relations  (4),  it  is  obtained: 


(7) 


— fP,)qr  +  c”(T-T#)+q,]  — 4»r*I  ^  =  0 


where  qr  is  given  by : 


qr  vp  MP 


2  v  RN£R  hR  2  Mp  lip 
R  P 


(9) 


(10) 


c.  Diffusion 

Assuming  that  only  concentration  diffusion  exists,  we  have : 


j  m; 


m 


4  n 


5e __L  ,  A  JL. 

r2l  \Y,  Yj  1J  Y,  dr  +  Y,  dr 


=  0  (11) 


where  «u  is  the  Lewis-Semenov  number,  given  by : 

°"= 

and  Y|  and  Djj  are,  respectively,  the  mass  fractions  and  diffusion  coefficients. 


(12) 


d.  Boundary  Conditions 

It  will  be  assumed  that  liquid  phase  reactions  do  not  exist  and  that  there  is  only 
one  liquid  chemical  species.  It  will  also  be  assumed  that  the  chemical  reaction  goes 
to  completion  at  infinity  and  that  the  temperature  and  composition  of  the  atmos¬ 
phere  at  infinity  are  known.  Thus,  we  have: 


r  =  r, 


T 

'i. 

fi» 


T. 

1 

0  (i  +  1) 


r  =  oo 


Too 

0  -*•  =  0 


(13) 


There  exists  one  boundary  condition  in  excess  of  the  number  cf  differential 
equations,  which  gives  the  burning  rate  m,  or  “eigenvalue”  of  the  system. 
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3.  SOLUTION  OF  THE  SYSTEM 

Solution  of  the  system  is  simplified  by  means  of  the  change  of  variables : 


6  =  -v-  ( 
Sr  \ 

T  — T„  +  J1-  ) 
cp  / 

mm  9 

(14) 

x  = 

cp  m 

AnX  r 

(15) 

From  which : 


d  e, 

dX 


fp 

7 


V  r,» 

X* 


wi 


de 

dX 


=-(0-*p) 


(16) 

(17) 


(lL- 

•i) 

1  a,,  4 — 

dV, 

1 

dYj 

[Mj 

lYi 

1  «  ^  y, 

dX 

dX 

where  X,,  given  by : 


(18) 


cp  jn 
AnX  r. 


(19) 


is  the  new  eigenvalue  of  the  system.  The  boundary  conditions  are  now  as  follows: 


X  =  X, 


e 

*i. 

'i. 


i 

0  (i*l) 


6  = 


(20) 


The  solution  of  the  problem  lies  on  the  integration  of  the  non-linear  system  of 
differential  equations  (16),  (17)  and  (18)  with  boundary  conditions  (20).  In 
addition,  an  expression  of  w,  as  a  function  of  the  mass  fractions  and  temperatures 
should  be  known. 

The  laws  of  variation  of  the  functions  *,  and  d  *,/dX  which  were  found  by  means 
of  numerical  integration  of  the  equations  for  several  typical  cases,  suggested  the 
adoption  of  an  approximate  analytical  integration  method.  This  method  is  based 
upon  considering  a  reaction  zone  of  finite  thickness  and  on  approaching  the  law  of 
variation  of  ^  within  such  zone  by  means  of  two  parabolic  curves  tangent  to  each 
other  at  point  (eloo  4-  eu)/2  and  tangent  to  the  lines  e{  =  *lao  and  e,  =  ejg  at  the 
boundaries  of  the  reaction  zone. 
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The  reaction  zone  is  determined  by  the  values  of  its  central  point  X*  and 
thickness  %.  The  maximum  value  of  d  f,/dX,  which  will  be  called  should  also 
be  fixed.  The  values  of  x,  c'*  and  X*  are  obtained  from  the  system  of  equations : 


(-*-) 


_  <P  X,2  r„2 


_  _p. 
I 


fi  =  f'i* 


X  =  X* 


z"|,  =  2|io,^dX  =  2<'“-''-)  (22) 

(w)  =0-X*w',*-4wi*  =  0  (2J) 

X  =  X* 

In  this  system  w  *  is  the  value  of  w,  (Y,,  8,  X)  for  X  =  X*,  in  which  Y,  and 
8  are  obtained  by  integrating  differential  equations  (17)  and  (18).  This  integra¬ 
tion  is  performed  by  introducing  into  these  equations  the  following  expressions 
for  f, : 

2  (-'~F---[X-(X«-  -|-»  J*  <X* - §-  <  x  <  X*) 


«l  =  V+2  (t|~  'u)  [X*  +  —  xjs  (X*  <  X  <  x*  +  i)  (25) 


In  regions  0  <  X  <  X*  —  */2  and  X*  4-  z/2  <  X  ^  X§,  e,  is  taken  constant 
and  equal  to  floo  and  respectively. 

The  integration  of  equation  (17)  is  straightforward.  X,  is  obtained  by  inte¬ 
grating  this  equation  between  ©t  and  8^.  The  integration  of  (18)  is  more  in¬ 
volved,  because  the  equation  is  not  linear.  However,  when  the  molecular  masses  of 
the  species  as  well  as  the  Lewis-Semenov  numbers  are  equal,  equation  (18)  becomes 
linear  and  the  integration  may  readily  be  performed. 

An  approximate  method  was  developed  to  to  integrate  Equations  (18)  for  the 
general  case.  This  method  is  based  on  the  series  expansion  of  the  expressions  for 
Yj,  assuming  that  x  >s  small,  which  is  the  case  that  normally  occurs  in  practice. 

Several  solutions  given  by  the  approximate  analytical  method  were  compared  to 
those  obtained  by  the  numerical  integration  of  the  equations.  Comparisons  were 
performed  for  a  wide  range  of  variation  of  the  characteristic  parameters  of  the 
process  and  for  different  reaction  rates.  These  comparisons  showed  that  the  approxi¬ 
mation  furnished  by  the  analytical  method  was  excellent  in  all  practical  cases. 


4.  THEORETICAL  RESULTS 

Chemical  kinetics  of  the  combustion  of  hydrocarbons  is  not  known.  Therefore,  in 
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order  to  perform  a  qualitative  study  on  the  combustion  of  fuel  droplets  in  air,  the 
combustion  process  is  approximated  by  means  of  a  second-order  over-all  reaction : 

A1  +  A2^2As  (26) 

and  chemical  kinetics  of  the  process  is,  then,  approximated  by  taking  for  the  reac¬ 
tion  rate  an  expression  of  the  form : 


Yi  Y2 
Mx  M2 


exp 


where  B  is  the  frequency  constant  and  E  the  activation  energy  of  the  over-all 
reaction. 

Expression  (27)  may  be  written  as  follows: 


^=-shfrQp!Y'Y2  <9-9»>-2“p( 

where  ©a,  ©0  and  Q  are  constants  given  by : 

F„E 


Rq, 

e0= 

0  q, 

C  2B 


Q  = 


Mr 


(28) 


(29) 

(30) 

(31) 


All  the  principal  variables  of  the  process :  burning  rate  m,  temperature  ©*,  flame 
radius  r*  and  maximum  temperature  ©max  are  obtained  as  a  function  of  the  dimen¬ 
sionless  product  yXp  r,,  being: 


A  = 


Qs 

MST 


(32) 


Burning  rates  are  generally  measured  by  the  slope  of  the  curves  of  the  square  of 
the  droplet  radius  as  a  function  of  time.  Such  nin  es  are  usually  approximated  by 
means  of  straight  lines,  and  the  slopes  are  proportional  to  X,,  because  they  are 
given  by : 


k  = 


m 

2  *  ei  rf 


(33) 


In  Figs.  2,  3  and  4  curves  of  X,  •  ©*,  ©max  and  r*/r§  are  represented  as  functions 
of  |/A  p  ra  for  several  values  of  Sx.  The  solutions  corresponding  to  the  case  of  an 
infinitely  fast  reaction  rate  are  also  included. 

All  results  lie  between  two  limiting  cases:  when  the  product  p  r(  tends  towards 
infinity,  the  values  for  Xg,  ©*,  ©max  and  r*/ra  tend  towards  those  obtained  by 
assuming  an  infinite  reaction  rate. 
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•  SURNINS  RATES  X, 


On  the  other  hand,  when  p  r,  tends  towards  zero,  the  value  of  X,  is  given  by : 

x.  =  'o«-^  04) 

r» 

9*  and  ©B1M  tend  towards  9^  and  r*/r,  tends  towards  infinity,  all  values  cor¬ 
responding  to  the  case  of  pure  evaporation  of  droplets  in  absence  of  combustion. 

Parameter  9%  influences  decisively  the  process.  When  ©,  is  small,  Xg,  6*,  0max 
and  r*/r§  are  close  to  their  asymptotic  values,  even  for  very  small  values  of  p  rB. 
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On  the  contrary,  for  large  values  of  such  variables  separate  from  their  asympto¬ 
tic  values  at  large  values  of  y7i  p  r,,  and  they  reach  rapidly  vaJues  very  close  to 
those  corresponding  to  y7i  p  r,  =  0.  (Ref.  10) . 

An  interesting  conclusion  is  that  solutions  close  to  the  case  of  pure  evaporation 
do  not  represent  actual  combustion  processes,  since  temperature  profiles  are  very 
flat  to,  with  maximum  temperatures  close  to  9 ec,  and  reaction  zones  are  very  wide. 
Therefore,  theoretical  results  show  that  minimum  values  of  either  the  pressure  or 
the  droplet  radius  exist  under  which  combustion  is  not  possible. 
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Variables  XB,  0*  and  r*/rM  change  rapidly  from  values  close  to  the  origin 
(j/X  p  rB  =  0)  up  to  values  close  to  those  corresponding  to  yX  p  r,  =  <».  Thefle- 
fore,  the  minimum  value  of  \Z~K  p  rB  required  for  combustion  may  be  approxi¬ 
mately  calculated  by  taking,  for  example,  the  value  of  ]/K  p  rg  corresponding  to 
the  average  value  of  X,  between  (X,)  _  Q  and  (X.)  ^  =  ^  .  Such 

minimum  values  depend  strongly  on  and  they  are  represented  in  Fig.  5. 

In  order  to  verify  the  aforementioned  theoretical  conclusions  an  experimental 
research  program  was  carried  out.  The  investigation  was  conducted  in  two  direc¬ 
tions:  in  the  first  place,  combustion  of  bromine  and  nitric  acid  droplets  in  gaseous 
hydrogen  was  observed.  Hydrogen  was  selected,  among  other  reasons,  because  the 
influence  of  chemical  kinetics  on  the  combustion  process  is  more  important  when 
the  molecular  mass  of  the  gaseous  species  surrounding  the  droplet  is  smaller  than 
the  molecular  mass  of  the  liquid  species  forming  the  droplet.  Results  of  this 
research  program  as  well  as  their  comparison  with  theoretical  results  have  already 
been  published  ih 12. 

Another  part  of  the  research  program  was  the  investigation  of  the  influence  of 
pressure  and  droplet  size  on  the  combustion  in  air  of  fuel  droplets. 

A  quantitative  agreement  between  experimental  and  theoretical  results  could  not 
be  expected,  due  to  the  lack  of  data  on  chemical  kinetics  of  the  combustion  of  fuels. 
However,  it  was  expected  that  at  low  values  of  the  product  p  r,  the  influence  of 
chemical  kinetics  should  decrease  the  slope  of  curves  r,2  =  f  (t)  and  it  should  also 
increase  the  values  of  the  ratio  of  the  flame  radius  to  the  droplet  radius.  Another 
important  objective  of  the  investigation  was  to  obtain  the  minimum  values  for 
combustion  of  the  product  p  rB. 

The  experimental  technique  consisted  in  observing  the  combustion  of  fuel 
droplets  suspended  of  thin  quarts  fibres.  The  droplet  size  could  only  be  modified 
within  a  narrow  range  (0.2 — 1  mm  in  radius).  Therefore,  the  influence  of  air 
pressure  was  the  main  objective  of  the  investigation.*) 

In  Fig.  6  the  research  facility  is  shown.  A  steel  combustion  chamber  was  used 
in  which  air  pressures  from  very  small  values  up  to  several  atmospheres  could  be 
established.  The  chamber  has  an  orifice  at  its  top,  which  is  closed  by  means  of  a 
glass  tube  which  moves  up  and  down  driven  by  an  electromagnetic  actuator.  When 
the  desired  pressure  is  set  within  the  chamber,  a  droplet  suspended  from  a  quartz 
fibre  is  introduced  into  the  tube,  closing  simultaneously  the  orifice.  The  tube  is 
then  moved  downwards  and  the  droplet  is  ignited  by  electrodes  which  are  placed 
in  position  by  the  motion  of  the  tube. 

Motion  pictures  of  the  burning  droplets  were  taken  by  means  of  an  Akeley 

*)  Hall  and  Diederichsen8  investigated  the  influence  of  pressure  on  the  combustion  of  fuel 
droplets.  However,  in  the  range  of  pressures  which  they  considered  (from  one  up  to  20 
atmospheres)  it  has  almost  no  influence  on  the  results,  because  for  sudi  conditions  results 
are  very  close  to  their  asymptotic  values.  Therefore,  they  concluded  that  pressure  did  not 
influence  appreciably  combustion  of  droplets. 
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camera,  and  the  values  of  both  droplet  and  flame  diameters  were  measured  as 
functions  of  time.  Droplet  diameter  was  taken  as  the  mean  value  of  two  perpendi¬ 
cular  diameters  measured  at  45°  from  the  fibre  axis.  Flame  diameter  was  taken  as 
the  value  of  the  maximum  horizontal  diameter  of  the  flame. 

Normal  heptane,  100  octane  gasoline,  RD  2498  kerosene  and  ethyl  alcohol  were 
tested.  Pressure  was  varied  for  each  fuel  from  the  minimum  value  required  to 
obtain  ignition  up  to  a  value  at  which  results  were  sensibly  constant  (independent 
of  pressure). 

Results  for  gasoline  and  kerosene  are  shown  in  Fig.  7  through  14.  It  may  be 
observed  that  the  experimental  laws  of  variation  of  the  evaporation  constants 
(proportional  to  X,)  and  of  the  flame/droplet  radius  ratio  r*/r,  as  functions  of  the 
pressure  and  as  functions  of  the  droplet  radius  agree  perfectly  well,  from  a  quali¬ 
tative  point  of  view,  with  those  obtained  theoretically.  It  may  be  pointed  out  the 
excellent  correlation  of  results  which  is  obtained  when  they  are  plotted  as  a  function 
of  the  product  p  r,  for  gasoline  and  as  a  function  of  p0-65  r,  for  kerosene. 

For  large  values  of  either  the  pressure  or  the  droplet  radius,  the  values  of  k  and 
r*/r,  tend  towards  asymptotic  or  constant  values,  and  when  the  pressure  or  the 
droplet  radius  are  reduced,  k  and  r*/r,  decrease  up  to  minimum  values  of  p  or  r, 
under  which  combustion  was  not  possible.  The  results  confirm  the  strong  influence 
of  chemical  kinetics  on  the  process  at  low  pressure  or  for  small  droplets  and  they 
imply  that  the  assumption  of  considering  an  infinitely  fast  reaction  rate  is  not  valid 
for  such  conditions. 
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Fig.  7 
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Fig.  8 


In  Figs.  15  and  16  photographs  taken  from  motion  pictures  of  burning  droplets 
of  gasoline  are  shown.  Fig.  15  shows  a  droplet  burning  at  ambient  pressure  and 
Fig.  16  a  droplet  burning  at  1/3  of  atmosphere.  In  this  case  the  flame  is  darker  and 
its  shape  is  almost  spherical,  which  means  that  free  convection  is  less  important. 
This,  in  turn,  implies  that  combustion  temperature  is  lower  according  to  the  results 
theoretically  predicted. 
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In  Fig.  17  and  18  minimum  values  of  the  droplet  radius  for  combustion  r80  are 
shown  as  a  function  of  the  pressure  for  the  four  fuels  investigated.  Such  radius  rto 
represents  the  minimum  value  for  which  ignition  could  not  be  achieved  by  means 
of  electric  sparks,  and  sometimes  the  extinction  value  of  the  radius  after  combus¬ 
tion,  of  the  type  shown  in  Fig.  19. 

It  may  be  expected  that  the  ignition  mechanism  should  influence  the  process  and 
there  also  is  a  possible  influence  of  the  flame  stability  on  the  extinction  values  of 
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Fig.  11 
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Fig.  12 

either  the  pressure  or  the  droplet  radius.  These  phenomena  explain  the  dispersion 
of  results  and  they  prevented  the  formulation  of  any  precise  law  of  variation  of  the 
extinction  radius  as  a  function  of  pressure.  However,  approximated  results  may  be 
obtained  from  Figs.  18  and  19. 

From  the  results  of  the  investigation  it  may  be  concluded  that  at  the  working 
conditions  which  exist  within  the  combustion  chambers  of  gas  turbines  or  compres¬ 
sion-ignition  engines,  the  droplets  of  their  sprays  will  be  under  conditions  at  which 
chemical  kinetics  will  influence  strongly  the  combustion  process.  Gas  pressure  is 
high,  but  the  droplets  are  small,  to  such  extent  that  the  product  pr„  will  probably 
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have  values  smaller  than  those  required  for  disregarding  chemical  kinetics.  It  is 
also  possible  that  for  certain  conditions  the  droplets  cannot  sustain  individual 
flames,  and  then,  the  combustion  process  would  consist  in  the  vaporization  of  the 
droplets  within  a  common  flame. 

Therefore,  spray  combustion  cannot  be  studied  by  applying  the  burning  laws  of 
individual  droplets  obtained  by  assuming  an  infinitely  fast  reaction  rate,  as  it  has 
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Fig.  16  C 
t  =  '*!*  sec. 


Fig.  16  B 
t  =  i/4  sec. 


Fig.  16  D 
t  =  alt  sec. 


Fig.  16  £ 
t  =  sec. 


Fig.  16  F 
t  =  ®/4  sec. 


Fig.  16.  A  gasoline  droplet  burning  at  1h  of  atmos. 
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been  customarily  done.  This  might  explain  some  of  the  discrepancies  between 
theoretical  and  experimental  results  which  usually  are  found  when  studying  spray 
combustion  and  it  also  gives  an  explanation  of  some  results  obtained  by  Graves  and 
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Fig.  19  A  tig.  19  B 

Ignition  Flame 


Fig.  19  C 
Extinction 

Fig.  19-  Combustion  of  a  n-hepttne  droplet  near  extinction  conditions  (p  =  XU  atmos.) 

Gcrstein  i3,  who  deduced  from  their  experiments  that  chemical  kinetics  should  in¬ 
fluence  strongly  combustion  of  sprays. 
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NOTATION 

A,  Q  Constants  of  the  reaction  rate  equation. 

B  Frequency  factor. 

Cp  Specific  heat  at  constant  pressure. 
cp  Average  value  of  cp  for  the  mixture. 

D[j  Diffusion  coefficient. 

E  Activation  energy. 

h|  Specific  enthalpy, 

k  Evaporation  constant. 

m  Burning  rate  (mass  of  fuel  burned  per  unit  time). 

Mj  Molecular  mass, 

p  Pressure. 

qi  Latent  heat  of  evaporation. 

qr  Heat  of  reaction, 

r  Radius, 

r,  Droplet  radius. 

R  Gas  constant, 

t  Time. 

T  Absolute  temperature, 

w,  Reaction  rate. 

X  =  Cpm/taAr,  dimensionless  coordinate. 

X,  =  cpl4jtlr,,  eigenvalue  of  the  system,  proportional  to  the  evaporation  constant. 
Y|  Mass  fraction. 

ej  Ratio  of  flux  of  mass  of  species  i  to  total  mass  flow. 

Thermal  conductivity. 

X  Average  value  of  X  for  the  mixture, 
ajj  =  XI Q  DjjCp,  Lewis-Semenov  number, 
vi  Stoichiometric  coefficients  of  chemical  reaction. 
q  Gas  density. 

Pi  Density  of  the  liquid  of  the  droplet. 

S  Dimensionless  temperature. 

©O  Constant. 

X  Reaction  zone  thickness  at  the  X  coordinate  system. 


SUBSCRIPTS 

i,  j  Chemical  species  in  general. 

R,  P  Reactant  and  products,  respectively. 

1  Gaseous  reactant  from  the  droplet. 

2  Oxidizer. 

3  Reaction  products. 

s  On  the  droplet  surface. 


SUPERSCRIPTS 

•  Center  of  the  reaction  zone.  (Point  at  which  d  fj/dX  is  a  maximum). 

.  Time  derivative. 

,  Derivative  in  respect  of  X. 

o  Standard  state. 
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